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Abstract: In a weak-coupling effective field theory framework we study quarkonium dis- 
sociation induced by inelastic scattering with partons in the medium. This is the dominant 
dissociation process for temperatures such that the Debye mass is larger than the binding 
energy. We evaluate the dissociation cross section and the corresponding thermal decay 
width. At leading order we derive a convolution formula relating the two, which is con- 
sistent with the optical theorem and QCD at finite temperature. Bound state effects are 
systematically included. They add contributions to the cross section and width that are 
beyond a quasi-free approximation, whose validity is critically reviewed. For temperatures 
such that the Debye mass is smaller than the binding energy, the dominant dissociation 
mechanism is gluo-dissociation consisting in quarkonium dissociation induced by the ab- 
sorbtion of a gluon from the medium. We calculate the gluo-dissociation cross section and 
width at next-to-leading-order accuracy. 
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1 Introduction 

Heavy- quarkonium suppression, first proposed in [1] as a signal of deconfinement, has been 
observed at SPS, RHIC and recently at LHC [2-5]. Although the current understanding 
is that the observed quarkonium suppression cannot be explained by cold nuclear matter 
effects alone [3] , the hot nuclear matter mechanism responsible for it is still under investi- 
gation. In [1], it was suggested that heavy quark-antiquark {QQ) bound states dissociate 
in a hot thermal bath because of colour screening of the QQ potential induced by the 
medium. In [6], another dissociation mechanism was identified in the Landau-damping 
phenomenon. Implications of the Landau- damping mechanism on the quarkonium dynam- 
ics, and anisotropic generalizations thereof, can be found in [7-10]. 

In the context of real-time thermal field theory, quarkonium in a thermal bath may be 
studied by taking advantage of the non-relativistic and thermal energy scales that charac- 
terize the system. The scales typical of a non-relativistic bound state are the heavy-quark 
mass m, the momentum transfer or inverse radius 1/r ~ mv and the binding energy 
E ~ mv'^. Because u <C 1 is the relative velocity of the heavy quarks, these scales are 
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hierarchically ordered: m 3> mv 3> mv'^. The thermal bath is charactetized by a tem- 
perature T and a Debye screening mass mo- In a weakly-coupled plasma, m/j ~ gT and 
also the thermal scales are ordered: T ^ niD- Integrating out systematically degrees of 
freedom associated with the highest energy scale leads to a hierarchy of low-energy ef- 
fective field theories (EFTs) [11, 12]. The ultimate of these EFTs can be interpreted as 
a finite-temperature version of potential non-relativistic QCD (pNRQCD) [13, 14]. Po- 
tential non-relativistic QCD describes the quarkonium dynamics through potentials and 
low-energy interactions. Thermal corrections affect both the real and the imaginary parts 
of the potentials. Thermal corrections to the real part of the colour-singlet potential may 
lead to the colour-screening phenomenon described by Matsui and Satz, whereas thermal 
corrections to the imaginary part of the colour-singlet potential induce a thermal width. 
Several mechanisms may contribute to the imaginary part of the potential. One of these is 
precisely the Landau-damping mechanism identified in [6]. In [11, 15], it was shown that 
the thermal width induced by the Landau-damping phenomenon is the principal source of 
quarkonium dissociation at weak coupling, responsible for keeping the quarkonium disso- 
ciated even at temperatures where colour screening on the real part of the potential has 
faded away. Moreover, lattice studies have shown that the potential may have a sizeable 
imaginary part also at strong coupling [16-18]. 

In a non-EFT framework, quarkonium decay widths induced by scattering with the 
medium constituents have been studied since long time (see for instance [19-22] and ref- 
erences therein). At leading order, two different dissociation mechanisms were identified: 
gluo- dissociation [23, 24] and dissociation by inelastic parton scattering [25, 26]. In the 
former case, the bound state absorbs a sufficiently energetic gluon of the medium and dis- 
sociates into an unbound colour-octet QQ pair. The gluon is physical, in the sense that 
its momentum is either light-like or time-like if it acquires an effective mass propagating 
through the medium. In the latter case a light parton of the medium, gluon or quark, 
scatters off the bound state by exchanging gluons, resulting again in its dissociation into 
an unbound colour octet. The momentum of the exchanged gluon is in this case space- like. 
In both cases, the decay widths were obtained by convoluting the T = cross section for 
the scattering process, possibly with some ad-hoc finite-temperature modifications, with 
the thermal distribution of the incoming parton. 

In the EFT framework, thermal decay widths have been investigated over a wide range 
of temperatures [11, 12, 27]. In [12], two mechanisms contributing at leading order to 
the quarkonium decay width were identified: singlet-to- octet thermal breakup and Landau 
damping. In the power-counting of the EFT, it was shown that the latter dominates over the 
former as long as m^) ^> E. These two mechanisms are not independent from those iden- 
tified earlier without EFT methods. In terms of elementary processes, the singlet-to-octet 
thermal breakup corresponds to gluon-dissociation and the Landau-damping mechanism to 
dissociation by inelastic parton scattering [28]. Beyond leading order the two mechanisms 
are intertwined and distinguishing between them becomes unphysical, whereas the physical 
quantity is the total width. 

The equivalence of the singlet-to-octet thermal breakup process to gluo-dissociation 
has been analyzed at leading order in [29]. There the singlet-to-octet thermal breakup 
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width and the gluo-dissociation cross section have been computed for several temperature 
regimes. The decay width in the regime mv ^ T ^ E ^ rriD, which has been suggested 
to be of relevance for T(IS') suppression at the LHC [30], agrees with the one previously 
calculated in [27]. Moreover it is consistent with the lattice QCD findings of [31]. The 
gluo-dissociation cross section agrees in the large- A'^c limit with the cross section derived 
in [32], which is often used in the literature. The gluo-dissociation cross section has been 
confirmed by [33] . 

In this paper, we will perform a similar analysis for what concerns the relation between 
the Landau-damping mechanism and the dissociation by inelastic parton scattering. We 
will compute in an EFT framework the dissociation cross section and the corresponding 
thermal width in different temperature regimes and relate the results with the existing 
literature. The cross section for dissociation by inelastic parton scattering was computed 
in [34] neglecting bound-state effects, and more recently (but in a different validity region) 
in [35]. The relevance of this process for heavy-ion collisions was realized in [25] and since 
then phenomenological expressions for the quarkonium decay width due to dissociation by 
inelastic parton scattering in the medium have been widely used in the literature (see [26, 
36-40]). As this process gives the dominant contribution to the decay width for rriD ^ 
E, which is the temperature regime at which quarkonium dissociates, validating those 
expressions from QCD is of utmost phenomenological relevance. In fact, we will show that 
the convolution formula commonly used in the literature to connect the dissociation cross 
section to the decay width does not follow from the optical theorem applied to QCD at 
finite temperature and we will suggest a different one. Furthermore, we will argue that 
neglecting bound-state effects is not a valid approximation in a weak-coupling setting below 
the quarkonium dissociation temperature. 

The paper is organized as follows. In section 2, we give some general arguments based 
on the optical theorem in thermal field theory and on effective field theories to derive 
a formula relating the dissociation cross section to the decay width. In section 3, we 
illustrate the basics of pNRQCD. In section 4, we study the cross section and width in the 
temperature regime T ^> mv ~ m^) and in section 5 we consider the regime T ~ mv ^ 
m^. In section 6, we focus on the temperature regime mv ^ T ^ m^ ^ E, where a 
multipole expansion in the temperature is possible. We obtain a colour dipole cross section 
and we study the validity region of this approximation. In section 7, we analyze the 
temperature regime mv ^ T ^ E ^ mn, where, as we shall show, both gluo-dissociation 
and dissociation by parton scattering get a contribution from the same EFT diagram. 
Hence, we compute next-to-leading-order (NLO) corrections to the gluo-dissociation cross 
section. Finally, in section 8 we draw some conclusions. 

2 General considerations on dissociation by inelastic parton scattering 

It has been suggested in [25] and used in most of the following literature on the subject 
that the width Fhq for dissociation by inelastic parton scattering of a heavy quarkonium. 
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HQ, at rest with respect to the thermal bath could be expressed by the convolution formula 



where the sum runs over the different species of incoming light partons, p, with momentum 
q = |q|, and the distribution functions, fp, are the Bose-Einstein distribution nB{q) = 
l/(exp(g'^/T) — 1) for gluons and the Fermi-Dirac distribution 'nF{q) = l/(exp(g^/T) + 1) 
for light quarks. The momentum qmin is the minimum incoming momentum necessary to 
dissociate the bound state. The distribution functions are convoluted with <t^^, a quantity 
identified with the parton-heavy-quarkonium dissociation cross section in the medium. 
In [25] and related literature, has been approximated by 2ap , where ap is the cross 
section of the zero-temperature process pQ — )• pQ. Because this approximation neglects 
bound-state effects, it is called quasi- free. ^ The scattering process pQ ^ pQ receives 
at leading order in perturbation theory contributions from the four diagrams shown in 
figure 1; these were computed in [34]. It is the cross section computed in [34] that is 
commonly used in eq. (2.1). Besides the distribution functions, additional thermal effects 
are usually added to (2.1) in the form of momentum-independent thermal masses affecting 
the dispersion relations and propagators of the light partons. Thermal masses also provide 
a cut off for the infrared divergences typically affecting the forward scattering amplitude. 



Figure 1. Tree level diagrams contributing to pQ —> pQ, where p is a parton and Q a heavy 
quark. Thick lines stand for heavy quarks, curly lines for gluons and thin lines for light quarks. 

Our purpose is to scrutinize eq. (2.1) and the quasi-free approximation in the light of 
QCD at finite temperature. We start by making some general considerations about the 
form of the dissociation formula. A quantitative treatment in a weak-coupling regime will 
be developed in the next sections. The dissociation process we are considering happens 
when partons in the thermal bath, i.e. partons with momentum and energy of the order 
of the temperature, scatter off the quarkonium and dissociate it by exchanging space-like 
gluons. First, we observe that by performing the calculation in Coulomb gauge, we only 
need to take into account transverse gluons, Ai, as external gauge fields. The reason is 
that in Coulomb gauge temporal gluons, Aq, do not thermalize at the scale T; moreover 
their spectral density vanishes and hence they do not contribute to cut diagrams. The 
Coulomb gauge is therefore a convenient gauge for the calculation, and we will adopt it 
for the rest of the paper. Next, we assume that the heavy-quark mass is much larger 
that the temperature: m ^ T. This implies that one can integrate out hard modes with 

^ In the literature, dissociation by inelastic parton scattering is sometimes called for conciseness quasi-free 
dissociation, after the approximation used for its computation. 
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energy and momentum of order m from QCD neglecting thermal effects, and replace QCD 
by non-relativistic QCD (NRQCD) [41, 42] as the fundamental theory. In NRQCD, the 
leading interaction between heavy quarks and gluon fields Ai is encoded in dimension-five 
operators; each of these interactions brings a suppression factor proportional to T/m. We 
conclude that in Coulomb gauge the last two diagrams in figure 1 are at least suppressed 
by a factor T/m with respect to the first two and can be neglected at leading order. 

Whenever the momentum carried by the gluon coupled to the heavy quark is larger 
than or of the same order as mv, based on the same argument given in the previous 
paragraph we can argue that the leading-order contribution of the first two diagrams in 
figure 1 comes from the coupling of a temporal gluon to the heavy quark. This case will be 
relevant for section 4 and 5. Whenever the momentum carried by the gluon coupled to the 
heavy quark is smaller than niv, but larger than the binding energy, the contribution of the 
temporal gluons cancels at leading-order in the multipole expansion in the square of the 
amplitude of the first two diagrams of figure 1 and the corresponding antiquark diagrams. 
This was noted in [12] and will be shown again in the following. Nevertheless, also in this 
case we may neglect at leading order the coupling of a transverse gluon to the heavy quark. 
The reason is that its contribution is proportional to the binding energy of the quarkonium, 
and this is subleading with respect to the momentum carried by the gluon. This case will 
be relevant for most of the temperature regions studied in the paper. Finally, whenever 
the momentum carried by the gluon coupled to the heavy quark is of the same order as 
the binding energy, the contribution coming from a transverse gluon coupled to the heavy 
quark will be at next-to-leading order in the multipole expansion as important as the one 
coming from a temporal gluon and we will need to consider them both. This case will be 
relevant only for section 7. 

Since in all cases only the first two diagrams in figure 1 contribute at non-vanishing 
leading order, it follows from the optical theorem that the dissociation cross section is 
proportional to the imaginary part of iD^'^ , where D^'^ is the gluon propagator. The 
gluon propagator may be computed in the so-called real-time formalism of thermal field 
theory (see e.g. [43]). A feature of this formalism is that the degrees of freedom double: 
while external particles are only of type "1", i.e. they live on the time-ordered branch of 
the Schwinger-Keldysh contour, in loops one has to consider also particles of type "2", 
i.e. particles located on the anti-time-ordered branch. It has been shown in [12, 44] that 
under the condition m ^ T heavy quarks do not thermalize up to exponentially suppressed 
contributions and can be treated as external probes. Hence, all vertices involving heavy 
quarks are of type "1" and, as a consequence, one needs to consider only the imaginary part 
of the "11" component of the gluon propagator coupled to the heavy quark. In general, this 
can be written in terms of an advanced (A), retarded {R) and symmetric (S) propagator 
as 

^fr(fco, k) = - [D'^^'iko, k) + D'X'iko, k) + D^,^{kQ, k)] , (2.2) 

where throughout the paper italic letters refer to the modulus of the spatial momentum, i.e. 

= |k|. Since (iDj^)* = iD^ , the discontinuity i-Dff — (z-Dj'J^)* is equal to iD'^'^ , which 
is in turn related to the retarded and advanced propagators through D^{kQ, k) = [1 + 
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2nB(|A;o|)] sgn(A;o) [Dj^{kQ, k) — {k^, k)]. The resummed retarded/advanced propagator 
is obtained by resumming the retarded/ advanced gluon self energy, n^*^^. The gluon 
polarization tensor satisfies relations similar to those valid for the gluon propagator. In 
particular, the "11" component of the self energy can be written as 

k) = ]^ [n^^(A;o, k) + n7(A;o, k) + (fco, k)] . (2.3) 

From (n^^)* = n^'^, it follows that the discontinuity 11 — (Ilj'j^)* is equal to H^^, since we 
also have that I\.^{ko,k) = [1 + 2nB(|A:o|)] sgn(A;o) {kQ,k) —Il^^{ko,ky\. In summary, 
the imaginary part of iD^'^ for spacelike momenta is proportional to iD^ , which, in turn, 
is proportional to n^^*^. In relation to the dissociation by inelastic parton scattering indeed 
only n^'' matters, because contributions coming from the poles of the propagator are not 
space-like. 

For the purpose of scrutinizing eq. (2.1), we will concentrate now on the case of a 
temporal gluon propagator in Coulomb gauge with incoming momentum k ^ kq. This 
is the only case needed in all sections of the paper with the exception of section 7. The 
resummed temporal gluon propagator in Coulomb gauge is given by 

It follows that the temporal symmetric gluon propagator can be written as 

which makes clear that the imaginary part of iD^ for space- like momenta comes from He?. 
At one loop, n^'^(A;o, k) is the sum of a light-quark, li^g ^{ko^ k), and a gluon contribution, 
g{ko, k). We will consider them for momenta k ^ k^. The symmetric gluon self 
energy can be computed by means of the cutting rules at finite temperature introduced 
in [45, 46] (see also [47, 48]). According to them, the quark contribution to the longitudinal 
polarization tensor reads 

Ufjk » fco) = Q dqq' (^1 - ^) nF(g)[l - nF{q)] , (2.6) 

where, here and in the following, we neglect contributions of order k'^/k'^ or smaller; is 
the number of light quarks in the self-energy loop. Noticing that —Tdnp{q) /dq = np{q) [1 — 
np{q))], eq. (2.6) can be rewritten as 

^fjk » ko) = ^^^-^ r dqqnpiq) , (2.7) 

which agrees with the expression of n^''^(^ 3> ko) that follows from eq. (45) of [12]. Si- 
milarly, the one-loop gluon contribution to the longitudinal polarization tensor in Coulomb 
gauge reads 

^s'gik » ko) = Q ^^'^^ " $ + ^) + ""^(^^l ' ^^-^^ 
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where Nc is the number of colours. Again, noticing that —TdnB{q)/dq = nB{q)[l + nB{q)], 
eq. (2.8) can be rewritten as 

" ^2 roo / 7,4 



^fjk » ko) = 



k^ f°° ( k^ \ 

-ns{k/2) + J^^^dqq^l-—^j nB(g) 



(2.9) 



which also agrees with the expression of Ilgg{k ^> ko) that follows from eq. (45) of [12]. 

In section 7 we will also need the symmetric transverse self-energy in the hard thermal 

loop (HTL) approximation ko,k <^T [49] . It reads 

j„2 / 1.2 \ /■°° 

n^(A:o, k) = ^eik'-kl) ( p - Ij j dqq^ iN,nB{q)[l + miq)] + Uf nF((z)[l - npiq)]) , 

(2.10) 

where 11"^ = {S^^ — k'^k^)W^ /2 and we have neglected higher-order corrections in k^/T. 

Equations (2.6), (2.8) and (2.10) suggest the following form for the parton-scattering 
dissociation width: 

ThQ = E / 7S3 U^) [1 ± ' (2-11) 

where the plus sign applies when the parton is a boson and the minus sign when the 
parton is a fermion. This expression incorporates the quantum-statistical effects both of 
Pauli blocking on the light-quark final states and of Bose enhancement on the gluon final 
states. We have set q^^ = q^^^ = q, which is a good approximation as long as T » 
since the incoming parton is on shell (hence qf^ = q) and its momentum is of the order 
of the temperature, while the transferred energy is of the order of the binding energy. 
Although eqs. (2.7) and (2.9) seem to allow for a dissociation width of the form (2.1), 
this is actually not the case if dp has to be understood as a cross section. In fact the 
quantity convoluted with the distribution functions, when using (2.7) and (2.9), cannot be 
interpreted as a parton-heavy-quarkonium cross section as it does not follow from applying 
the optical theorem to a quarkonium-quarkonium amplitude. Our conclusion is therefore 
that eq. (2.1) in its common interpretation is not justified by QCD at finite temperature. 
QCD at finite temperature suggests instead formula (2.11) or its generalization for the case 
9in 7^ lout- the rest of the paper, we will explicitly derive eq. (2.11) at leading order for 
a wide range of temperatures. 

The momentum qmin is equal to the absolute value of the quarkonium binding energy. 
As long asT ^ E, which will be the case for all thermal regimes discussed in this paper, we 
can set Qmin = in the convolution integral (2.11). Corrections in qrain/T are suppressed. 

As a final comment, we remark that eq. (2.11) only holds in a leading-order picture, 
where one light parton with momentum of order T scatters off the bound state. At higher 
order, when more partons appear in the initial or final states, or when some of them have 
momenta of order gT, eq. (2.11) is no longer valid. 

3 General considerations on the EFT approach 

Before calculating the dissociation width and cross section in an EFT framework, we sum- 
marize here some general aspects of this framework. As mentioned in the introduction. 
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this approach is based on the hierarchies of non-relativistic and thermal scales typical 
of quarkonium in a quark- gluon plasma. For definiteness, we will assume that the non- 
relativistic scales, the temperature and the Debye mass are larger than the typical hadronic 
scale Aqcd; which justifies a perturbative treatment for all of them; this also implies that 
V ^ Og- A system whose energy scales may possibly satisfy this assumption is the bot- 
tomonium ground state at LHC [30]. 

In the following, we will consider several possible temperature regimes at weak cou- 
pling. We will proceed integrating out from thermal QCD all scales larger than E ~ mv'^. 
The ultimate EFT that describes QQ pairs with momentum of order mv and energy of 
order mv'^ interacting with gluons and light quarks of energy and momentum of order mv'^ 
or smaller has the form of pNRQCD. If the temperature is larger than mv'^, the matching 
coefficients of pNRQCD will depend on the temperature. At the accuracy we will need it 
in the following, the pNRQCD Lagrangian at weak coupling reads [13, 14] 

>CpNRQCD = Cught + y TV {st [ido -hs]S + Qt [iDo - K] O 

+0^r ■gBS + S^r-g'EO + ^O^ {r ■ gB , O}} , (3.1) 

where £iight is the part of the Lagrangian that describes the propagation of light quarks 
and gluons, S = S U/VK and O = \/20" T" are the QQ colour -singlet and colour-octet 
fields respectively, E is the chromoelectric field and iDoO = idoO — gAoO + OgAo- The 
trace is over colour and spin indices. Gluon fields depend only on the centre-of-mass 
coordinate and on time. In (3.1) we have neglected irrelevant operators of order r^, 1/m 
or smaller; we have also neglected quantum corrections to the matching coefficients of 
the dipole operators, which are of order or smaller, and beyond our accuracy. In the 
centre-of-mass frame, the singlet and octet Hamiltonians have the form (p = — iVr): 

h,^^=P-+Vj'J + ^ + \ + ..., (3.2) 

m ' m 

where the dots stand for higher-order terms in the 1/m expansion. The first two terms 
in the right-hand side, which are the kinetic energy and the static potential respectively, 
constitute the leading-order Hamiltonian. 

We refer to [12, 44] for details on pNRQCD in the context of the real-time formalism 
of thermal field theory. We only mention that the "2" QQ fields decouple, hence all singlet 
and octet QQ fields appearing in pNRQCD amplitudes have to be understood as "1" fields. 



4 The T ^ mv ^ mD case 

In the same framework that we adopt here, the case T mv ^ mj:) was studied in 
section VI of [12] and in section V B of [11]. Since m ^ T, we start by integrating 
out from QCD the mass scale, obtaining NRQCD. Then we integrate out the scale T 
from NRQCD to arrive at a version of NRQCD that is modified by the temperature. In 
particular, the gauge and light-quark degrees of freedom are described by the hard thermal 
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loop Lagrangian [50, 51], whereas the heavy-quark sector is not modified at leading order. ^ 
The next step consists in integrating out also the scales mv and m^) to obtain a version of 
pNRQCD specific for the hierarchy T ^ mv ~ m£).^ This version of pNRQCD corresponds 
to the Lagrangian (3.1), with Cng^t given by the HTL Lagrangian and with thermally 
modified QQ potentials. In particular, vj'^^ is at leading order the potential computed 
in [6]; it reads 

(r) = -C,aJ'— + mn) + ^2C, a. T T dt ( '^^^^^ - l) 
\ r J Jo \ mort J 

where Cp = {N^ - l)/(2iVc) and 



{t^ + 1)2 ' 

(4.1) 



-l, = ^fiVc + ¥V (4.2) 



2 J 

The imaginary part of eq. (4.1) describes precisely the physics of dissociation by inelastic 
parton scattering. Corrections of higher order in 1/m and g are beyond the accuracy of 
this paper. 

At leading order in the multipole expansion, the equation of motion for the singlet 
field resulting from the pNRQCD Lagrangian is a Schrodinger equation with the potential, 
provided by eq. (4.1), consisting of a real and an imaginary part. As discussed in [12], if 
mv ^ mo the latter is larger than the former by a factor of T/mo and the bound state 
can be considered dissociated. On the other hand, if mv is sufficiently larger than mo, 
corresponding to the situation T ^ mv ^ m/5, then both the imaginary part and the 
screening are perturbations of the Coulomb potential [11, 12]. The temperature at which 
the real and imaginary parts become of the same size can be defined as the dissociation 
temperature T^. One then has ~ mg^/^ [11, 15] (see also [53] for numerical estimates of 
the T(IS') dissociation temperature). 

We derive now the cross section and decay width from the potential (4.1) under the 
assumption that the real part of the potential is larger than its imaginary part. At leading 
order, the decay width is given by 



nl 



-(n,/|21myj°)(r)|n,/) , (4.3) 



where |n, /) is an eigenstate of p^/m + 'R.eVs^\r), and n, / are the principal and orbital 
angular momentum quantum numbers identifying the quarkonium. From eq. (4.1), it 
follows that 



TT Jo {t^ + V tr J 



(4.4) 



^ This thermal version of NRQCD was called NRQCDhtl in [44, 52]. In the abelian case the corre- 
sponding EFT was named NRQEDt in [53] . 

^ This version of pNRQCD was called pNRQCDhtl in [52] . In the abelian case the corresponding EFT 
was named pNRQEDt in [53]. In [44] it was instead named pNRQCDmD- We also remark that in the 
literature there exist two versions of pNRQCDhtl, one for T mv, corresponding to the present case, 
and one for mv S> T, which was studied in detail in [27, 54]. 
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We observe that, given eq. (4.2), the imaginary part of the potential can be separated 
in a part coming from the scattering with light quarks and in a part coming from the 
scattering with gluons. More in detail, following the arguments of section 2, the imaginary 
part originates from the symmetric part of the longitudinal propagator taken in the HTL 
limit, kQ,k ^ T. The longitudinal propagator is given by eq. (2.5). The longitudinal 
polarization tensor, n|gP(A;o, k), in the HTL limit follows from eqs. (2.6) and (2.8) expanded 
for ko,k <C T. For the quark contribution we have 

- ^nO%(A;o = 0, « T) = J j^riFiq) [1 " Md)] = "^Dlquark , (4.5) 
where m|)|quark = g^nfT"^ /Q. For the gluon contribution we have 

- 2^n00^(^0 = 0, « T) = ^ j ^nB(g) [1 + n^{q)] = ml\^,^,^ , (4.6) 
where mj^lgiuon = Q^NcT"^ /2>. We can then rewrite the width as 

nF(g)(l - nF(g))(n,/|Sg(r,g)|n,/) + nB{q){l + nB{q)){n,l\T.g{r,q)\n,l) , 

(4.7) 
where 

and 

. , ^ o ("^ dtt ( sin(tr)\ , , 
Finally, we can identify 

<(g) = (n,/|Sp(r,g)|n,0, (4.10) 

with the cross section of a quarkonium state (with quantum numbers n, /) with a parton 
p = q, g in the medium, and arrive at the formula (2.11). Note that the gluon- and quark- 
induced cross sections differ only in the colour structure. 

We conclude this section with some comments about the cross sections. First, we note 
that T,p{q,r) and hence the cross sections (7p\q) do not depend on the parton momentum, 
q ^ T. This holds at leading order for momenta such that m ^ q ^ mv. As can be seen 
from eqs. (2.6) and (2.8), it is a consequence of assuming the temperature to be much larger 
than the other scales, mu, rriD and E. We also remark that the cross sections depend on the 
temperature only through the Debye mass, rriD- Hence, in contrast to what happens for the 
gluo-dissociation cross section [29], we cannot relate the cross section for inelastic parton 
scattering, not even at leading order, with a zero temperature process. The underlying 
reason is the infrared sensitivity of the cross section at the momentum scale mv, which can 
be seen by putting m/j = in eqs. (4.8) and (4.9). This infrared sensitivity is cured by the 
HTL resummation at the scale mo, as it will become more apparent in the next section. 
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4.1 The mv ^ mu case 



We consider now the cross section and decay width of a quarkonium state in the special case 
mv ^ mi). Under this condition the state is Coulombic, i.e. KeVj'^\r) ~ —Cpag/r. This 
case may be possibly realized only for a quarkonium IS" state [30], which we are going to con- 
sider in the following. The IS* wave function is given by (rllS*) = l/(-y/7raQ^^) exp(— r/ao), 
where oq = 2/ [mCpOis) is the Bohr radius. The evaluation of the dissociation cross section 
follows then easily from 

^^^|Sin(tr),_, 16 



tr ' ' {f^al + Af 
In particular, the light-quark cross section becomes 



= (15'|Sg(r,g)|lS') = SttCfUj a^alfimDao 



where 



and a very similar formula holds for the gluon cross section: 



|l-g)= u,.2 , (4-11) 

(4.12) 

2 r - 16 + 8x2ln(4/3;2)" 
/(x) = -^ 1-4 ^ , (4.13) 

(4.14) 
(4.15) 

acg = d,-KCFNcalal, (4.16) 

so that CTp'^ = (JcpfimDao), with p = g,q. Plugging eqs. (4.12) and (4.14) into eq. (4.7), 
we obtain the width, which reads 



erf (g) = {lS\^g{r,q)\lS) = STrCpN^ataifimDao) 



It is convenient to define the constants 



acq = STrCpUf ttg , 



and 



r,s = 2CFa,T 



^{mnaoY - 16 + 8(mDao)^ln (4/(771000)^) 
((m£,ao)2 - 4)3 



(4.17) 



The width in this regime was already obtained in eq. (1.7) of [55]. Our result appears to 
be larger by a factor of 2 due to the fact that in [55] the width is defined as one half of 
ours. Under the assumption mv ^ mo we can further expand f{m£,aQ) for m^aQ <C 1, 
obtaining for the cross section 



and for the decay width 



cp 



In 



m%al 



CpdsT mj^QQ ( In 



m%al 



with p = q, g, 



(4.18) 



(4.19) 
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Figure 2. Cut diagrams in NRQCD. The momentum of the cut partons is of order T. Dashed hnes 
stand for the cuts, thick hnes with arrows for heavy quarks and antiquarks, curly hnes for gluons 
and grey blobs for either light-quark or gluon loops. 



I 




Figure 3. Cut diagram in pNRQCD. The momentum of the cut partons is of order T. The single 
line stands for a QQ pair in a colour-singlet configuration, the double line for a QQ pair in a colour- 
octet configuration, the circle with a cross for a chromoelectric dipole vertex, the curly line for a 
HTL gluon and the grey blob for either a light-quark or a gluon loop. 

An alternative derivation of the cross section and decay width given in eqs. (4.18) and 
(4.19) consists in assuming right from the start that T ^ mv ^ and in evaluating 
the potential through the hierarchy of EFTs introduced in section VI E of [12]. It was 
shown there that the potential, as defined in the ultimate EFT, receives a contribution 
from the scale mv and one from the scale rriD- The leading-order real part of the static 
potential is the Coulomb potential, which comes from the scale mv. For what concerns 
the imaginary part of the static potential, the contribution from the scale mv is infrared 
divergent and originates from the cut diagrams in figure 2 when the momentum flowing 
in the gluon is of order mv. It reads in dimensional regularization (D is the number of 
spacetime dimensions) 

Imy™^(r) = ^a,r''Tml (^^^ + 7£; + In^ + ln(r ^)2 - 1^ , (4.20) 

where je is the Euler-Mascheroni constant and /i is the subtraction scale. The contribution 
from the scale mo arises from the cut diagram in figure 3, where the displayed gluon stands 
for a HTL-resummed gluon carrying a momentum of order m^- The gluon interacts with 
the QQ field through the chromoelectric dipole interactions induced by the (second line 
of the) Lagrangian (3.1). The diagram was evaluated in eq. (87) of [12] and found to be 
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ultraviolet divergent. It gives 

luiV^^ir) = -^a,r^Tml(-^-jE + ln7r + ln^ + l) , (4.21) 

which holds at leading order in an expansion in E/ttid.'^ At this order the octet propagator 
in figure 3 can be taken to be l/(— /cq + ie), where is the gluon momentum; this means 
that the rescattering of the unbound final-state heavy quarks can be neglected. Summing 
the two contributions the divergence and related fj, dependence cancel and one obtains 




Figure 4. The rescaled cross sections (Jp^ /fJcp as a function of mrtCLQ- The continuous line corre- 
sponds to firriDao), as defined in eq. (4.13), whereas the dashed line follows from eq. (4.18) and is 
thus the expansion of f{moao) for m/^oo <C 1. 

In figure 4 we plot the rescaled cross section jocp^ obtained in eqs. (4.12) and 
(4.14), as well as its expansion for itidclo 1, obtained in eq. (4.18). They are the 
continuous and dashed curve respectively. The two curves overlap up to mzjao ~ 0.2. We 

* The condition T S> mv implies nio ^ E for a Coulombic bound state in a weakly-coupled plasma. 
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note that the curves have physical meaning only for values of m£,aQ significantly smaller 
than 1. The reasons are that the dashed curve has been obtained as the leading term of an 
expansion in itidclq, and the continuous curve, although following from the imaginary part 
of the potential (4.1), which is valid also for moaQ ~ 1, assumes the bound state to be 
Coulombic, see eq. (4.11), which is instead a valid assumption only for m£,aQ much smaller 
than 1. 

4.2 Comparison with the literature 

As we have already mentioned, dissociation by inelastic parton scattering was first consid- 
ered in the context of heavy-ion collisions in [25]. There the cross section was approximated 
by twice the parton heavy-quark scattering cross section computed in [34], an approx- 
imation called quasi-free approximation, complemented by the addition of momentum- 
independent thermal masses. The only dependence on the bound-state properties comes 
through gmin, which was obtained from an in- medium binding energy calculated from a 
phenomenological potential model. 

The parton heavy-quark scattering cross section corresponds to the square of the first 
two diagrams in figure 1. Conversely, the optical theorem relates the cross sections that we 
have computed in section 4 with the square of the sum of the diagrams shown in figure 5. 
Each of these diagrams individually corresponds to those computed in [34], however, once 
the square is taken, interference terms appear in the form of diagrams with gluons attached 
to diff^erent heavy-quark lines. These terms are the ones that are sensitive to the spatial 
separation between the heavy quark and antiquark, and ultimately to the bound state. 
The quasi- free approximation consists in neglecting these interference terms. 



Figure 5. Diagrams contributing to the (amplitude of the) dissociation cross sections derived in 
section 4 and given by eq. (4.10). The thick lines represent heavy quarks and antiquarks while the 
thinner lines represent light quarks from the medium. 

More precisely, in the right-hand side of eq. (4.4), the term proportional to sin(tr) /{tr) 
is the interference term while the remaining term comes from the last two cutting diagrams 
shown in figure 2. These two distinct terms may be traced in the cross sections. For 
instance, in the Coulombic case, mv ^ rriD^ we have that, in the right-hand side of 
eq. (4.13), the term is the contribution from twice the square of the parton heavy- 
quark scattering diagram, whereas the remaining terms, which are neglected in the quasi- 
free approximation, are the interference terms and give a non-trivial dependence on the 
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wave function and thus on the properties of the bound state. In figure 6 we plot for 
comparison jacy = /(m^joo) and the quasi-free term 2/ {niDaQY . For small values of 
mijao, where (4.12) and (4.14) are valid, the quasi-free approximation overestimates the 
cross section by more than one order of magnitude. Conversely, for large values of m/jao, 
where the quasi-free approximation is sensible, the two curves overlap. In this region, 
however, we cannot treat the quarkonium as a Coulombic bound state, and, therefore, 
eqs. (4.12) and (4.14) are no longer valid. 




0.1 0.2 0.5 1.0 2.0 5.0 10.0 



Figure 6. The rescaled cross section a^^ /acp as a function of mnao- The continuous line corre- 
sponds to f{mDO-Q), as defined in eq. (4.13), whereas the red dashed line corresponds to 2/{muao)^, 
which is the only term that survives in the quasi-free approximation. 

It is furthermore important to note that under the condition mv ^ niD the parton 
heavy-quark scattering contribution to the dissociation cross section is exactly cancelled at 
leading order in m/^ao by a contribution coming from the interference terms. The remnant 
of this cancellation, which is entirely due to interference/bound-state terms, is what makes 
up for eq. (4.18). The cancellation was noticed in [12, 56]. 

In summary, the quasi- free approximation is justified when all relevant thermal scales 
and, in particular, are much larger than mv. However, in a weak-coupling regime, 
the condition mo ^ mv requires a temperature T ^ mg, which is parametrically larger 
than the dissociation temperature T^^ ~ mg^^^. On the other hand, whenever m^ 
mv, not only the quasi-free approximation is not justified, but its contribution is exactly 
cancelled by bound-state effects. We conclude, therefore, that, at least in a weak-coupling 
framework, the quasi-free approximation is not justified for all range of temperatures where 
a quarkonium can exist. 
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5 The T ~ mv ^ mn case 



The temperature region T ~ mv ^ niD was studied in detail for the case of muonic 
hydrogen in [53]. In the muonic hydrogen case, which may be interpreted as an abehan 
version of heavy quarkonium, thermal corrections are encoded in electron loops. In [53] 
both the real and the imaginary part of the potential were computed. Here we focus on 
the imaginary part of the potential, which is the quantity relevant for dissociation; thermal 
corrections are encoded in light-quark and gluon loops. 

Let us briefly sketch our procedure. Since m ^ T we can use NRQCD as a starting 
point. Then, as T ~ mu, these two scales have to be integrated out at the same time. 
In doing so we go from NRQCD to a new particular thermal version of pNRQCD. In 
the light-quark and gluon sector it is made of the HTL Lagrangian. In the heavy-quark 
sector the Lagrangian is as in (3.1) with the potential of a form specific to the hierarchy 
T ~ mv ^ m£). It is at leading order a Coulomb potential that receives small corrections 
with a complicated functional dependence on the temperature. A feature of these thermal 
corrections is that they are infrared divergent. The infrared divergence that we observe 
in the imaginary part of the potential cancels against ultraviolet divergent contributions 
coming from the scale mn, exactly as in the previous section. We discuss first light-quark 
and then gluon contributions. 

The leading light-quark contribution to the imaginary part of the static potential comes 
from the one- loop fermion contribution to the diagrams in figure 2. At this order only the 
coupling of the heavy quarks with temporal gluons is relevant. Hence, we may follow the 
discussion of section 2, and the contribution reads 

^ 4 V -) -'-^"^^ 

A f Q / \ r / \n A D f ^ k G ^ f ^ \ / \ 

= 2vr5 C,nf j j^.n^m " n^{<l)] ^ ' J^^^^ (2vr)^-i I' " V J ' ^'"'^ 

where we have used the symmetric longitudinal polarization tensor given in eq. (2.6). The 
integral in k has an infrared divergence that we have regulated in dimensional regular- 
ization. This divergence cancels against contributions coming from the scale m^- The 
contributions from the scale can be evaluated from the cut diagram in figure 3 when 
the momentum flowing in the loop is of the order of the Debye mass. Since E ~ ma^ is 
parametrically smaller than mn (exactly by one power of g), we can expand the intermedi- 
ate octet propagator for E -^mo neglecting rescattering effects and the result is the same 
as in eq. (4.21). Because that equation sums the contribution of quarks and gluons into 
m\), we need to use eqs. (4.2) and (4.5) to disentangle the one from the other. Summing 
the quark contribution from the scales T ~ mv with the quark contribution from the scale 
mD the divergences cancel and the light-quark contribution to the imaginary part of the 
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potential reads 



144 V3 



+ 



(2^ 



nF(g)[l - nF(g)] 



/ sin(ir) 



2q 



V tr 



4tg^ 



tr 



(5.2) 



The decay width follows from eq. (4.3). It can be written in the form of eq. (4.7), with 



Sg(r,, 



Svr 



4 , , ^ 6 dt / sm( tr) 

3 /"^q /sin(tr) ^ 



(5.3) 




Figure 7. Comparison between the rescaled cross sections a]^ jucq for T ^ mv ^ mu and 
T ~ rnu ^ m-D, corresponding to eqs. (4.12) and (5.4) respectively. The dashed red line is the 
function /(mj^ao), the black continuous line is the function hq(mDao, 10) and the dotted blue line 
is hqimoaQ, 1). 

For T ~ mv ^ uid, the cross section of a weakly-coupled 15 quarkonium state with 
light quarks from the medium is given by 



f^l^il) = crcghq{mDao,qao 



where 



hq{x,y) 



In 



2 \l + y^ 



2^1n(l + .^), 



(5.4) 
(5.5) 
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Gcq is defined as in eq. (4.15) and we have made use of the matrix elements computed 
in eqs. (4.11), (4.23) and (4.24). In the previous section we derived that, in the case 
T ^> mv ~ rriD, crg^{q) = crcqf{m£)ao), with / given in eq. (4.13). For T ~ mv ^> ttid, 
and differently from the case T ^ mv ^ m£,, the cross section depends on the momentum. 
However, since T ~ mv ^ m£) imphes that m^ao <^ 1 and qaQ ~ 1, in the hmit mDO-o — s- 
and qaQ — t- oo the functions / and hq should coincide. This can be seen in figure 7, where we 
have plotted hq as a function of m/jao for qao = 10 (for larger values of qao the plot would 
not change significantly). The residual difference between /(mijao) and hq{mDao, 10) for 
moderate values of mnao is due to the fact that in computing / we have resummed the 
HTL interaction while in computing hq we have not; the difference is however a small 
perturbation as long as moo-Q ^ 1- 

The leading gluon contribution to the imaginary part of the static potential can be 
computed similarly to the quark contribution. From the scale T, we have 



ImF/ (r) 



It 4- 
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(5.6) 



where we have used the symmetric longitudinal polarization tensor given in eq. (2.8). As 
in the quark case, the integral in k has an infrared divergence that cancels against the 
contribution coming from the scale m£). The contribution from the scale m/j is given 
in eq. (4.21), whose gluonic part can be disentangled by means of eqs. (4.2) and (4.6). 
Summing the gluon contribution from the scales T ~ mv with the gluon contribution from 
the scale mjy the complete leading-order gluon contribution to the imaginary part of the 
potential reads 



ImK 
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d'^q 



TnB{q)[l + riBiq)] 
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The decay width follows from eq. (4.3). It can be written in the form of eq. (4.7), with 
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(5.8) 



/o \ " / 2g^r^ \^ (gr)2 

For T ~ mv ^ mpi, the cross section of a weakly-coupled 15* quarkonium state with 
gluons from the medium is given by 



15 



(q) = (rcghg{mDao,qao) , 



(5.9) 



-18- 




0.0 0.1 0.2 0.3 0.4 0.5 

Figure 8. Comparison between the rescaled cross sections cr]f /<Jcg for T ^ mv ~ rriD and 
T ^ mv :s> mo, corresponding to eqs. (4.14) and (5.9) respectively. The dashed red hne is the 
function /(m/jao), the black continuous line is the function hglmuao, 10) and the dotted blue line 
is hgimoao, !)■ 

where 

h,(..y) = - 1„ (^) - 5 + + 1.. ( J^) - ^ Ml + v') . (5.10) 

and acg is defined as in eq. (4.16). To obtain eq. (5.9), besides eqs. (4.11), (4.23) and (4.24), 
we have made use of the matrix element 

(isf-^ns) = (5.11) 

In figure 8 we plot hg and / as a function of muao, and show that for large values of 
qao and low values of niDO-o the cross sections calculated in this section and in section 4 
overlap. 

The imaginary parts of the static potential given in eqs. (5.2) and (5.7) have been 
calculated here for the first time. We observe that the imaginary part of the potential at 
the scale T, given by eqs. (5.1) and (5.6), incorporates both the effects of the free-quark- 
parton scattering in the r-independent part and of the bound state in the r-dependent 
part. In contrast, the contribution to the imaginary part of the potential coming from the 
scale niD, given by eq. (4.21), is just an r-dependent contribution, hence it would be set 
to zero in the quasi-free approximation. 

6 The mi; > T > m/) > E case 

The temperature region mv ^ T ^ mu ^ E was studied in the static case in [12]. Since 
thermal contributions associated with higher-order terms in the 1/m expansion turn out 
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to be suppressed, at leading order we may indeed restrict ourselves to the static case. Our 
procedure goes as follows. Because mv ^ T we can use the pNRQCD Lagrangian at T = 
as our starting point. Next we integrate out the scale T to define a version of pNRQCD 
specific to the hierarchy mv ^ T (see footnote 3). Its Lagrangian was derived in [27, 54] 
and it features thermal corrections to the potentials in the heavy-quark sector and the 
HTL Lagrangian as /^ught- The correction to the colour-singlet QQ static potential has 
an infrared divergence that is compensated, as in the cases previously discussed, by an 
ultraviolet divergence coming from the scale niz). After integrating out mu, we obtain a 
new version of pNRQCD whose potential incorporates corrections coming from the Debye- 
mass scale. The new potential is finite and renormalization-scheme independent. All the 
effects discussed in this section are specific of having an interacting QQ system and would 
be absent in the quasi-free approximation. 

The effect of the light-quark loop at the scale T comes from the diagram in figure 3 
when the momentum flowing in the loop is of order T. Only the longitudinal part of the 
gluon propagator contributes and yields [12] 
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where we have used that n5j'(A;o, k) +Yi^^{kQ, k) is real and even in fco- The octet propagator 
appears in (6.1) as l/(— fco + ie), which is again a consequence of working at leading order 
in E/T <^ 1 and neglecting rescattering effects. It is because the even part of l/(— feo + ie) 
is proportional to 6{kQ) that only longitudinal gluons contribute, for the chromoelectric 
dipole interaction due to transverse gluons is proportional to kQ. Equation (6.1) would 
also follow from expanding eq. (5.1) in r. As in that case we have regulated the infrared 
divergence in dimensional regularization. The contribution from the scale m^) is the one 
computed in eq. (4.21), whose quark and gluon contributions may be disentangled by means 
of eqs. (4.2), (4.5) and (4.6). The reason is again that E/mo ^ 1 and we are working at 
leading order in E/mo- The fact that we are working at leading order in E/T and E/rriD 
is ultimately also the reason why 1/m effects provide only subleading thermal corrections 
to the static result. Summing the contributions coming from the scale T and the scale 
mu we get a finite expression for the light-quark contribution to the imaginary part of the 
potential, which, cast in (4.3), provides a thermal decay width of the form (4.7) with 
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From this it follows that the cross section of a weakly- coupled 15" quarkonium state with 
light quarks from the medium reads 
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Figure 9. Rescaled dissociation cross sections due to light quarks, a]^ jocq-, as a function of the 
rescaled momentum qa^. The dashed blue curve displays the cross section for mv 3> T ^ rriD 3> E 
given in eq. (6.3), the continuous black curve displays the cross section for T ~ mv ^ niD given 
in eq. (5.4), and the dot-dashed red curve displays the cross section for T ^ mv ^ mjj given in 
eq. (4.18). For all the curves we have assumed mnao — 0.1. 



where acq has been defined in eq. (4.15). 

Performing a similar calculation for the part involving the gluon loop, at the scale T 
we obtain 
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The contribution from the scale mjj is the gluonic part of eq. (4.21); summing it to eq. (6.4) 
we get the gluonic contribution to the imaginary part of the potential, which, cast in (4.3), 
provides a thermal decay width of the form (4.7) with 
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From this it follows that the cross section of a weakly- coupled 15" quarkonium state with 
gluons from the medium reads 

4q^- 



where acg has been defined in eq. (4.16). 
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Figure 10. Rescaled dissociation cross sections due to gluons, a)^^ juc-g-, as a function of the rescaled 
momentum qa^. The dashed blue curve displays the cross section for mv ^ T ^ mo ^ E given in 
eq. (6.6), the continuous black curve displays the cross section for T ~ mv ^ mo given in eq. (5.9), 
and the dot-dashed red curve displays the cross section for T ^ mv ^ m^ given in eq. (4.18). For 
all the curves we have assumed mnao = 0.1. 

We can now compare the results obtained in section 5 for T ~ mv m/5 to the two 
limiting cases: T ^ mv 3> mo, discussed in section 4.1, and mv ^ T ^ m^ » E, 
discussed here. Physically this means going from temperatures close to the dissociation 
temperature, T^, to temperatures in which all the interactions with the medium can be 
described at leading order by a chromoelectric dipole vertex. In figure 9, we plot the three 
cross sections with light quarks for m^ao = 0.1. The continuous black curve is the result 
for T ~ mv ^ m^ as given by eq. (5.4): for qao ^ 1 it is indeed approximated very well 
by the dashed blue line, which is the result for mv ^ T ^ m^, ^ E just obtained in 
eq. (6.3); for qao > 4 the black line is well approximated by the constant, dot-dashed red 
line, which is the result for T ^ mv ^ mo, as given in eq. (4.18). In figure 10, similar 
curves show the three cross sections with gluons. 

6.1 Comparison with the literature 

The results of this section bear a direct relation to the momentum diffusion coefficient k 
of a single heavy quark, first introduced and computed in [57]. According to the field- 
theoretical definition of [58], the diffusion coefficient can be written as the time integral of 
the thermal expectation value of two chromoelectric fields linked by Wilson lines stretching 
along the temporal axis. Our eqs. (6.1) and (6.4) correspond indeed to the integral of the 
correlator of two chromoelectric fields up to a factor r"^. Equation (B 13) of [57] matches 
the structure of eqs. (6.1) and (6.4), while eq. (B 14) of [57] contains the same g-dependent 
factor, [ln(4g-^/m|)) — 2], that we find in eqs. (6.3) and (6.6). Only in this section, when 
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the thermal scales are set between the bound state scales mv and mv'^, which clearly do 
not appear in the single quark case, is this direct comparison possible. It would then be 
also possible to use the NLO computation of k in [59, 60] to obtain some of the order g 
corrections to the width, while bound-state dependent contributions of possibly the same 
order, like those coming from the expansion of the octet potential, would need a new 
dedicated computation.^ 

7 The mt; > T > > case 

The temperature region mv ^ T ^ E ^ mo was studied in detail in [27]. This tem- 
perature region is technically more difficult to treat than the other ones discussed in the 
paper. One of the reasons is that we find diagrams whose imaginary parts contribute both 
to gluo-dissociation and dissociation by inelastic parton scattering. Another reason is that 
the calculation of these diagrams involves not only longitudinal gluons but also transverse 
gluons. Since the temperature is smaller than the typical inverse radius of the quarkonium, 
all quarkonium interactions with the medium are described at leading order by chromoelec- 
tric dipole vertices. Hence, also the effects discussed in this section are specific of having 
an interacting QQ system and would be absent in the quasi-free approximation. 

Our starting point is the pNRQCD Lagrangian defined in the previous section after 
integrating out the temperature T. The reason we can start from there is that the La- 
grangian of an EFT is only sensitive to the hierarchy of energy scales above its ultraviolet 
cutoff. Therefore integrating out the scales mv and T is unaffected by the relation between 
mD and E. Because the Lagrangian at the scale T is the same, the leading contributions 
to the imaginary part of the static potential from the scale T can be read off directly from 
eqs. (6.1) and (6.4). 

The next step consists in integrating over momenta of the order of the binding energy E. 
At this scale, the octet propagator in the diagram of figure 3 can no longer be taken as 
l/{—k() + ie), for the momentum /c^ flowing in the loop is of the same order as the octet 
energy. This means that the rescattering between the unbound heavy-quarks and their 
relative motion has to be taken into account. As a consequence, the interaction between 
heavy quarks and transverse gluons, which is proportional to /cq, does not vanish and 
contributes to the computation. This should be contrasted with what occurs for other 
energy scales. 

The retarded /advanced propagator of a transverse gluon in Coulomb gauge after HTL 
resummation reads 



^ We remark that the NLO calculation of [59, 60], when applied to heavy quarkonium, shows how the 
distinction between dissociation by inelastic parton scattering and gluo-dissociation fails beyond leading 
order. For instance, the intricate structure of cuts discussed in [59, 60] includes processes with two light 
partons in the initial and final state. 
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where Il^^{ko/k) = {5'^^ — yy)W^{kQ ± ie,k)/2 is the HTL retarded/ advanced trans- 
verse gluon self energy.^ It has the following properties. First, n]J^(0) = 0. Second, 
Imll^ ^(/co/A;) is different from only for |A;o//c| < 1. The fact that the imaginary part 
of ^ does not vanish only for space-like momenta implies that it contributes only 
to quarkonium dissociation through inelastic parton scattering. Finally, we have that 
ReH^^ (jto/k > 1^ ~ rn'jj. The equation k^ — k"^ — Rell^J^ {k^/k) = has then a solu- 
tion only for time-like momenta. This solution is called plasmon pole and obeys a special 
dispersion relation with a momentum-dependent mass. For time-like momenta the imag- 
inary part of ^ comes only from the ie prescription on the plasmon pole. Hence the 
plasmon pole contributes only to the quarkonium gluo-dissociation. There are two distinct 
momentum regions where the gluon propagator is nearly singular or singular. One is the 
region where ko ^ k ^ E and /cg — A;^ ~ E'^. In this region, which has been called off-shell 
region in [27] , we can expand the gluon propagator in the transverse gluon self energy. The 
other is the region where ko ^ k ^ E and A;g — /c^ ~ m|). In this region, which has been 
called collinear region in [27], HTL effects have to be resummed. 

The situation at the scale E is therefore the following. Gluons interact with the QQ 
pair at leading order through chromoelectric dipole interactions. The relevant Feynman 
diagram is shown in figure 3. Because at this scale the rescattering of the QQ pair cannot be 
neglected, both longitudinal and transverse gluons contribute to the quarkonium thermal 
decay. Quarkonium may decay either by emitting time-like or light-like gluons, which 
corresponds to cutting the diagram in figure 3 along the gluon propagator and picking up 
its pole contribution, or by scattering with partons in the medium. This last situation 
corresponds to cutting the diagram in figure 3 along the gluon self-energy diagram and 
picking up its discontinuity, which is encoded in the symmetric polarization tensor. The 
momentum of the gluon is in this case space- like. According to our definitions the first decay 
process contributes to quarkonium gluo-dissociation whereas the second one to dissociation 
by inelastic parton scattering. Both decay processes are intertwined at the scale E and may 
be disentangled only by looking at the time-like or space- like nature of the gluon interacting 
with the QQ pair. If the gluon interacting with the QQ pair is longitudinal, then the 
residue of its plasmon pole contribution is exponentially suppressed for momenta /cq ~ ^ ~ 
E ^ mo [63, 64], whereas a contribution to the thermal width comes from the imaginary 
part of the longitudinal polarization tensor. This is different from zero only for space- 
like momenta and hence contributes only to quarkonium dissociation by inelastic parton 
scattering. If the gluon interacting with the QQ pair is transverse, then it may contribute 
either through its plasmon pole or through the imaginary part of the transverse gluon 
self energy. The former case, which may happen only in the collinear momentum region 
for time-like gluon momenta, contributes to quarkonium gluo-dissociation. The latter case, 
which may happen both in the collinear and in the off-shell momentum regions for space- like 

^ The explicit expression of is [61, 62] 

feo /fco ,\feo, f ko + k±ie \' 
fc2 J 2k '^\ko-k±ieJ ■ 



n?,A(fco/fc) 



rriD 
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momenta, contributes to quarkonium dissociation by inelastic parton scattering. We will 
disentangle the gluo- and parton-scattering dissociation contributions to the cross section 
in the following two sections. 

7.1 Gluo-dissociation 

Quarkonium gluo-dissociation was studied in [29] at leading order in an uid/E expansion, 
which corresponds to evaluating the diagram in figure 3 with a free gluon. Here we add 
HTL effects, which amounts at computing the gluo-dissociation width and cross section at 
NLO. It is when HTL effects are taken into account that also parton-scattering dissociation 
happens. 

For the reasons discussed in the previous paragraphs, gluo-dissociation is due to the 
diagram shown in figure 3 when the gluon is transverse and its momentum time-like. The 
relevant contributions have been calculated in appendix A of [27]. Using those results the 
gluo-dissociation thermal width at next-to-leading order in mo/E and E/T reads 

Tni = ^asCpTin, l\ri{E - /if )Vi|n, /) + '^a.Cpin, l\ri{E - /if )Vi|n, /) 



^ -{fi^^n 



In 



l{E-h^^Y 
ml 



ri\n,l), (7.2) 



where h^^ = p^/m + a^/ {2Ncr) is the leading-order octet Hamiltonian. The first term 
in the right-hand side of (7.2) is the leading (zeroth-) order width in the rriD/E and E/T 
expansions: it reproduces the result of [29] for T ^ E. The other two terms are the 
next-to- leading-order corrections . 

The gluo-dissociation width may be expressed as the convolution of a gluon-heavy- 
quarkonium dissociation cross section in the medium, cr^l^^, and a gluon distribution func- 
tion: 

'^''^ m{q)a^,M- (7.3) 



- nl 



Note that, in contrast to the parton-scattering dissociation case, described by eq. (2.11), 
there is just one parton of the medium involved in the gluo-dissociation process and there- 
fore just one distribution function appearing in (7.3). Comparing (7.3) with (7.2) and 
expanding the Bose-Einstein distribution for T ^ E, we obtain the gluo-dissociation cross 
section 

a^UQ) = Z{q/mn)4^;^{q), (7.4) 

where cr^l^^^ is the leading-order cross section that corresponds to the first line in the right- 
hand side of (7.2). Its explicit expression for a 15 Coulombic bound state can be found 
in [29, 33] and reads 

4p 



= ^2'».V(P + 2)=^ (t(,f + p=) . (7.5) 



where p = 1/{N^ — 1), t{q) = ^yq/\El\ — 1 and Ei = —mC'^al/A is the energy of the 
first Bohr level. The absolute value of Ei provides the low-momentum cut-off in the 
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Figure 11. IS* gluo-dissociation cross section in units of ao/m. The continuous black line is 
the leading-order result, o'gj^^''^ The dashed blue line shows o-gj^^, according to eq. (7.4), for 
mo mag = 0.2 and the dotted red line for rriD mag = 0.5. 

integral (7.3). The factor Z can be understood as a wave- function normalization of the 
gluon due to the HTL resummation, it reads 

Z{x) = 1 - 1^ [ln(8x2) - 2] . (7.6) 

The effect of the normalization factor Z on the IS gluo-dissociation cross section is 
shown in figure 11, where the cross section is expressed in units of ao/m = 2/ {m^CFCts) 
and the gluon momentum in units of \Ei\. The plot shows how the HTL resummation 
results in a global lowering of the cross section.^ 

7.2 Dissociation by inelastic parton scattering 

Contributions from the scale E to quarkonium dissociation by inelastic parton scattering 
come from the different sources that we have analyzed in the introduction of section 7. 
We have the contribution from longitudinal gluons and that from transverse gluons. The 
contribution from transverse gluons is divided into contributions from the collinear region 
and from the off-shell region, which are separated by a cut-off. Only the sum of all these 
contributions is gauge invariant and cut-off independent. These different contributions have 
all been computed: the contribution from the longitudinal gluons can be found in eq. (5.17) 

^ Z{q/mD) becomes larger than one for q < 2''^^^emD ~ 0.96mD. However the cross section has a 
threshold at q = \E-i\, and | _Ei | is larger than mo in the assumed hierarchy. 
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of [27]; the contributions from the transverse gluons, of which we have to keep only the 
contributions coming from space-hke momenta, can be found in appendix A of [27]. 
The final hght-quark loop contribution to the decay width from the scale E is 



nl, q 



37r 



(n, /|r^|n, /) 



rnF(g)(l - nY{q)) 



1 
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+ 



IE 
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+ In 



ln(27r) 
J 



(7.7) 



Adding to it the contribution from the scale T, as given by eqs. (6.1) and (4.3), the 
divergence cancels and we can cast the decay width in the form (4.7) with 



2 ^2 



IGnCFUf r 
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For a IS Coulombic state, the corresponding quark-heavy-quarkonium dissociation cross 
section then reads 



cq 



Inl ^ 1 + ln2-2 



m 
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(7.9) 



where (Jcq has been defined in (4.15). 

The gluon loop contribution to the decay width from the scale E is the same up to a 
different colour structure and different distribution functions: 
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Adding to it the contribution from the scale T, as given by eqs. (6.4) and (4.3), the 
divergence cancels and we can cast the decay width in the form (4.7) with 
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For a 15" Coulombic state, the corresponding gluon-heavy-quarkonium dissociation cross 
section then reads 

.2 ■ 
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+ ln2 - 2 



(7.12) 



where Ucg has been defined in (4.16). 

The parton-scattering decay width is of order a^T x {mD/mvY , therefore suppressed 
by a factor (mz)/i?)^ with respect to the gluo-dissociation width, which at leading order 
scales like agT x {E/mv)'^ . The parton-scattering decay width is comparable in size to the 
next-to-leading-order correction to the gluo-dissociation width that appears in the second 
line of (7.2), while the next-to-leading-order correction appearing in the first line of (7.2) 
is of order a^T x [E/mv)'^ x (E/T). Hence, in the temperature region mv ^ T ^ E ^ 
ruD, dissociation by inelastic parton scattering is a subleading effect with respect to gluo- 
dissociation and may be neglected in first approximation. At the same time we observe 
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that in all other temperature regions examined in the paper we had mu ^ E. Therefore, 
in those regions, just the opposite holds and dissociation by inelastic parton scattering is 
the parametrically dominant quarkonium dissociation process. These observations agree 
with the early findings of [12]. The dominance of dissociation by inelastic parton scattering 
over gluo-dissociation at high temperatures was noticed in [25, 26]. 

7.3 Comparison with the Uterature 

In [37] parton-scattering and gluo-dissociation have been treated in an unified framework. 
We will highlight some qualitative features of that work that are common also to other 
approaches but that are different from the EFT treatment presented here. The first dif- 
ference is that the calculation of [37] uses the formula (2.1) for both the parton-scattering 
and the gluo-dissociation widths. We have seen that this is consistent with QCD only in 
the latter case. The cross sections used in [37] have been derived from a Bethe-Salpeter 
framework in [35]. The calculation includes systematically bound-state effects, but with 
some limitations: it is valid in the large Nc limit, hence it neglects rescattering effects of the 
unbound colour-octet quark-antiquark pair; it describes the quarkonium interaction with 
gluons through chromoelectric dipole vertices, hence the description holds for gluons whose 
energy and momentum are smaller than mv. The cross section does not include systemati- 
cally thermal effects, for it is calculated at zero temperature. Constant thermal masses have 
been added to regulate infrared divergences. This amounts at a phenomenological tuning: 
from a QCD perspective one should recall that momentum and temperature-independent 
masses are neither consistent with HTL resummation nor with weak-coupling perturbative 
calculations. 

8 Conclusions 

Quarkonium dissociation through scattering with light partons is one of the processes 
responsible for the thermal decay width of quarkonium in a deconfined medium. It is 
the dominant process for temperatures such that the Debye mass, m^i, is larger than the 
binding energy, E. 

We have studied this process in a weak-coupling effective field theory framework, where 
quarkonium dissociation through scattering with partons in the medium may be related, for 
momentum transfer larger than E, to the imaginary part of the potential and the Landau 
damping phenomenon. We have shown that in our setting the quasi-free approximation, 
which consists in approximating the dissociation cross section of the quarkonium with that 
of two free quarks, is never a valid approximation. In particular, for momentum transfer 
smaller than or of the same order as the inverse radius of the quarkonium, the dissociation 
cross section in the quasi-free approximation is exactly cancelled by bound-state effects. 

The parton-scattering dissociation cross section, valid for temperatures such that rriD 
is much larger than E, is of the form (4.10), with Sg given in (5.3) for scattering with 
light quarks from the medium and given in (5.8) for scattering with gluons. In the 
specific case of a Coulombic 15 state, the cross sections with gluons and quarks are given 
respectively by eqs. (5.4) and (5.9) and shown by the black curves in figures 9 and 10. In the 



-28- 



region of temperatures where mo is of the order of the inverse radius of the quarkonium, 
i.e. in the region where screening effects are important, T,q is given by eq. (4.8) and Tig 
by eq. (4.9). In the case of a Coulombic 15 state, the cross sections with quarks and 
gluons are given respectively by eqs. (4.12) and (4.14). The impact of the screening on 
the cross sections is shown in figures 7 and 8. The parton-scattering dissociation cross 
section has been computed also for temperatures such that m^) is much smaller than E. 
The light-quark contribution is given in eq. (7.8) with the corresponding cross section for 
a IS Coulombic state in (7.9), and the gluon contribution is given in eq. (7.11) with the 
corresponding cross section for a 15 Coulombic state given in (7.12). In this temperature 
regime, dissociation by inelastic parton scattering is, however, a subleading effect. 

From the parton-scattering dissociation cross section we may calculate the correspond- 
ing dissociation width through eq. (2.11). This expression is justified by general arguments 
based on the optical theorem and cutting rules at finite temperature, and by the explicit 
calculations performed in the paper in the different temperature regimes. Equation (2.11) 
should replace the widely used formula (2.1), which is justified only in the gluo-dissociation 
case. 

Gluo-dissociation is the process occuring when a sufficiently energetic gluon of the 
medium is absorbed by the quarkonium and dissociates it into an unbound colour-octet 
QQ pair. This is the dominant dissociation process in the temperature region where mu is 
much smaller than E, while it is subleading with respect to dissociation by inelastic parton 
scattering if m/j is much larger than E. Gluo-dissociation has been called singlet-to-octet 
thermal breakup in the effective field theory literature on the subject. We have calculated 
the gluo-dissociation cross section in eq. (7.4) and thermal width in eq (7.2) at next-to- 
leading order in itid/E and E/T. This is currently the most accurate determination of 
gluo-dissociation in weak coupling, whose impact is shown in figure 11. 
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